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First order superconducting phase transition in CeCoIns. 
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We investigated the magnetic field dependence of the superconducting phase transition in heavy 
fermion CeCoIns (T c = 2.3 K) using specific heat, magneto-caloric effect, and thermal expansion 
measurements. The superconducting transition becomes first order when the magnetic field is ori- 
ented along the 001 crystallographic direction with a strength greater that 4.7 T, and transition 
temperature below To ~ 0.31T C . The change from second order at lower fields is reflected in strong 
sharpening of both specific heat and thermal expansion anomalies associated with the phase tran- 
sition, a strong magnetocaloric effect, and a step-like change in the sample volume. The first order 
superconducting phase transition in CeCoIns is caused by Pauli limiting in type-II superconductors, 
and was predicted theoretically in the mid 1960s. We do not see evidence for the inhomogeneous 
Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) superconducting state (predicted by an alternative theory 
also dating back to mid-60's) in CeCoIns with field H \\ 001. 

PACS number(s) 74.70. Tx, 71.27.+a, 74.25. Fy, 75.40. Cx 



The behavior of superconductors in magnetic field has 
primary scientific and technological importance. It un- 
derlies such diverse areas as magnetic imaging, energy 
transmission and storage, ultrasensitive instrumentation 
and electronics, and many other fields of technology and 
medicine. At the same time it reflects a very fundamental 
property of matter - the behavior of electrons in magnetic 
field. BCS theory, presented in 1957 (Ref .111) , gave micro- 
scopic explanation of a number of phenomena observed 
during the previous half century of research on supercon- 
ductivity. The theories put forth in the early 1960's, that 
addressed the effect of magnetic field on superconductiv- 
ity, were the first extensions of BCS that made predic- 
tions of new phenomena and provided tests of BCS the- 
ory's predictive powersoTj. Magnetic field can suppress 
superconductivity via two effects: orbital pair breaking of 
superconducting pairs in the superconducting state and 
Pauli paramagnetism due to electron spins, which low- 
ers the relative energy of the normal state. It was shown 
that when the Pauli effect is sufficiently strong relative to 
the orbital effect, the superconducting phase transition 
may change fjjPjm second order (BCS result for zero field) 
to first ordeioti This is due to competition between two 
energies basic to the understanding of superconducting 
and normal states of metals: condensation energy of su- 
perconducting pairs and magnetic energy of the normal 
electron spins due to Pauli paramagnetism. This predic- 
tion is very straightforward, and yet eluded confirmation 
for almost 40 years. This is partly the reason why the ex- 
perimental observation of the first order superconducting 
transition became one of the "holy grails" of the super- 
conducting community. 

A number of conventional superconductors were pro- 
posed as candidates for observation of the first order su- 
perconducting transition in magnetic field, due to their 



high orbital critical field H C 2Q (weak orbital pair break- 
ing) and, therefore, relatively strong Pauli limiting ef- 
fect, in the early and mid-sixties. Kpaserimental search, 
however, came up with null resultsoi!3. The failure to 
observe the first order superconducting transition was 
attributed to high spin-orbit scattering rate in all of the 
compounds investigated^. Here we present specific heat, 
magnetocaloric, thermal expansion, and magnetostric- 
tion data for CeCoIns that demonstrate that the su- 
perconducting phase transition indeed changes from sec- 
ond to first order at a critical point with temperature 
To = 0.31T C , in very good agreement with theoretical es- 
timate. Recently, a first order phase transition at low 
temperature in CeCoIns was inferred from the step in 
thermal conductivity in CeCoIns at H C 2 || [001], which 
was suggested to be "likely due to an entropy jump"Ell. 
Our specific heat data offer a direct and unambiguous 
proof of the first order phase transition in CeCoIns with 
magnetic field in this orientation. 

CeCoIns is a recently diseavered ambient pressure 
heavy fermion superconductorO with the record high su- 
perconducting transition temperature T c = 2.3 K for this 
class of compounds. A number of thermodynamic and 
spectroscopic measurements indicate spin-singlet, even 
pairing state, path lines of nodes in the superconducting 
energy gapE^TtJ. Therefore, the lowest possible symme- 
try of the superconducting state of CeCoIns is d-wave. 
Recently, Izawa et alO reported the four-fold modulation 
of thermal conductivity of CeCoIns in magnetic field in 
support of the d x 2_ y 2 order parameter, similar to high 
temperature superconductors. 

First order superconducting phase transition in 
CeCoIns occurs in magnetic field close to the supercon- 
ducting critical field H c2 = H c2 {T = 0) = 4.95 T, with 
field along the [001] crystallographic direction. The 
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FIG. 1. Specific heat and entropy of CeCoIns. (a) Spe- 
cific heat vs. temperature. Closed symbols - decay method, 
numbers over the peaks indicate magnetic field in T. Open 
symbols - heat the pulse method. (A) 4.62 T, (y) 4.72 T, 
(o) 4.77 T, (<) 4.8 T, (o) 4.87 T, (>) 4.925 T. (b) Entropy as 
a function of temperature from (a). From left to right: 4.925 
T, 4.8 T, 4.75 T, 4.7 T, 4.6 T, 4.5 T. Arrows indicate step-like 
features in entropy at T c for H > 4.7 T. 

change from the second order nature of the transition, 
observed at zero and low magnetic field, to first order at 
high field occurs at T « 0.72 K = 0.31T C . We used a 
variety of techniques to investigate the region of interest 
of the H — T phase diagram, including measurements of 
specific heat at constant magnetic field, magnetocaloric 
effect, thermal expansion, magnetostriction, and resistiv- 
ity. 

Fig. |l|(a) shows specific heat as a function of tempera- 
ture for several values of magnetic field in the high field 
region. The data represented with open symbols were 
collected with the standard heat pulse technique. When 
the phase transition becomes first order, it sharpens sub- 
stantially. In this regime the temperature decay method, 
where specific heat is extracted directly from the tem- 
perature trace of the system coming to equilibrium, was 
particularly useful to resolve the specific heat anomaly 
(solid symbols). The two data sets for H = 4.8 T, ob- 
tained with the temperature decaying down (stars) and 
up (dots), as well as the data collected with heat pulse 



method (left open triangles) overlap each other, indicat- 
ing good internal thermal equilibrium of the cell. 

The specific heat anomaly for the second order phase 
transition (e.g. 4.5 T data) displays a characteristic step 
at T c (predicted to be AC/C = 1.43 in the BCS theory), 
and then gradually drops below T c . In CeCoIns AC/C 
dron^ steadily from the very high value of 4.5 at zero 
fieldB to w 1.1 at H = 4.5 T. The temperature width 
of the anomaly at half maximum AThm scaled by T c for 
4.5 T is u;(4.5 T) = AT H m/T c | 4 .5T » 25%. 

As we move into the regime of the first order phase 
transition, the maximum of the specific heat anomaly 
rises, e.g. AC/C(4.8T) = 1.9, and it becomes substan- 
tially narrower, w(4.8T) = 6%. It should be kept in 
mind that this sharpening of the anomaly takes place in 
the region of the phase diagram where the boundary be- 
tween the normal and superconducting states is crossed 
at more of a glancing angle during the temperature sweep 
as H — ► H C 2- As a result, the anomaly should broaden 
if it remains second order, in contrast to the evolution of 
the data. 

Figure |l|(b) shows the entropy S of CeCoIns as a func- 
tion of temperature for various values of the magnetic 
field, obtained from the data in Fig. |l|(a). There is a 
clear difference in the temperature evolution of the en- 
tropy for the field below and above Hq ss 4.7 T. Below 
this field (T c > 0.7 K) transition manifest itself by a kink 
in S. For fields above 4.7 T there is a step-like feature in 
S at T Cl as expected for a first order phase transition. 

To elucidate the low temperature behavior and deter- 
mine the temperature To of the critical point (and corre- 
sponding field Hq) at which the superconducting tran- 
sition changes from second to first order, with higher 
precision, we studied the phase diagram of CeCoIns by 
sweeping magnetic field. This was done under close to 
adiabatic (constant entropy) conditions, regulating the 
bath temperature to be approximately that of the sam- 
ple, which was weakly thermally coupled to the bath. 
The behavior of the system is then governed by the mag- 
netocaloric effect. 

Figure ||(a) shows several sweeps of magnetic field up 
and down, starting at different temperatures. Below H 
there is a sharp change in temperature as magnetic field 
crosses the phase boundary. Temperature drops when the 
phase boundary is crossed as the field is swept up, since 
the system goes from the low entropy to the high entropy 
phase, and the temperature has to decrease to keep the 
entropy constant!!!]. The temperature swing is reversed 
(temperature rises) on the down field sweep. The change 
in the temperature of the sample after the crossing of the 
first order phase boundary AT p b is a measure of the latent 
heat associated with transition. T c corresponds to the 
maximum of the derivative ^ , and AT p b is determined 
by extrapolations of the fits to the H vs. T curves outside 
of the transition region to T c . An example of such a 
procedure is displayed in the inset of Fig. |2|(a) for the data 
with T c = 0.41 K, where horizontal dotted line segment 
represents AT p b- The difference AT p b should equal zero 
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FIG. 2. Quasi-adiabatic sweeps of magnetic field, (a) H 
vs. T during magnetic field sweep for various starting (T, H) 
points. Dotted (dashed) line is a phase boundary in second 
(first) order transition region. Arrow - same as in (b). In- 
set: magnetocaloric field sweep with T c = 0.41 K < To and 
T c = 0.74 K > To. Dotted lines are cubic fits to the data 
and AT p i, for data with T c — 0.74 K. (b) Change in temper- 
ature at the phase transition during magnetic sweeps in (a), 
(o) field swept up, (A) field swept down. Arrow indicates 
the temperature To at which the superconducting transition 
changes from first to second order. Dotted line - guide to the 
eye in the region of the first order transition close to To. 

for the second order phase transition. This is indeed ob- 
served, within experimental scatter, for the data with 
T c > To, illustrated in the inset of Fig. ^(a) for the data 
with T c = 0.74 K. Figure ||(b) displays the measured 
ATpb as a function of T c . The cross-over from the first 
order transition with non-zero AT p b to the second order 
transition with AT p f, = occurs at a sharply defined crit- 
ical temperature To = 0.72 ± 0.02 K, which is indicated 
by the arrow. 

We also observed the change from the second to 
first order nature of the superconducting transition of 
CeCoIns via thermal expansion measurements, depicted 
in Fig g. The coefficient of the thermal expansion, 
a{T) = l~ 1 dl/dT, was determined down to 50 mK by uti- 



lizing an ultrahigh-resolution capacitive dilatometer with 
a maximum sensitivity of Al/l > 10 . Thermal expan- 
sion along the crystallographic [001] direction was mea- 
sured for two different plate-like single crystals in mag- 
netic fields up to 8 T applied along [001]. For one of the 
crystals isothermal magnetostriction measurements were 
performed as well at T = 0.2 K and at T = 1.5 K. 

In the low field range H < 4 T a step-like anomaly 
in a, indicative of a second-order transition, is observed, 
which shifts towards lower temperatures upon increasing 
H. With increasing fields (see Fig. ||) the signature in 
a sharpens anomalously and becomes peak-like, with ex- 
tremely high absolute values of a. This again indicates 
a change of the nature of the superconducting transition 
from second order to first order for magnetic field on the 
order of 4.6 T. A first order transition should result in 
a jump in the sample length, corresponding to a diver- 
gence of a. Thus the peak-like signature in a indicates 
a broadened first order transition. Additional evidence 
for this interpretation comes from the isothermal magne- 
tostriction experiments displayed in the inset of Fig. 3. 
Whereas at T = 1.5 K the kink in Al/l at H = 3.55 T 
indicates a second-order phase transition, the jump in 
Al/l, observed for T = 0.2 K at H = 4.86 T, provides 
clear evidence for the first-order nature of the transition. 

How do our experimental results compare with the- 
oretical predictions? Pauli paramagnetism leads to an 
upper limit for the magnetic field— H p = Ao/v2/iB) 
called Clogston paramagnetic limittH, which supercon- 
ductor can support. Here A is the superconducting 
energy gap, and \ib is the Bohr magneton. Orbital 
effects of magnetic field also limit JT C 2. The relative 
strength of the orbital pair breaking of magnetic field 
and Pauli limiting can be characterized bv. the parame- 
ter a = \/2H C 2o/Hp, introduced by MakiB, where JT C 20 
is an orbital critical field m the absence of the Pauli lim- 
iting. Maki's calculations!] show that for a > 1 the sec- 
ond order phase transition between the normal state and 
Abrikosov vortex state becomes unstable and changes to 
first order at a higher field, and if orbital effect is ne- 
glected (a — oo), this change takes place at the reduced 
temperature t = T /T c = 0.55.0 

Also in the early 60's, an alternative theory, also based 
on Zeeman energy of electron spins in magnetic field, 
suggested that a new spatially inhomogeneous super- 
conducting state, now called the Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state, may be stabilized close to the 
criticaLfield H c % at which superconductivity is suppressed 
to zeroM. FFLO state forms a wedge between the nor- 
mal state and the homogeneous Abrikosov vortex state. 
The tricritical point where the phase boundaries between 
the three states meet denotes the point of instability of 
second order phase transition at low field to the appear- 
ance of the FFLO wedge above the critical rfield within 
the FFLO scenario. Gruenberg and Gunthem (GG) gen- 
eralized the FFLO theory and included the orbital effect 
of magnetic field. The authors conclude that when the 
orbital pairbreaking effect is sufficiently small (a > 1.81), 
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3. Thermal expansion coefficient a vs. T of CeCoIns 
for fields 4T < H < 4.9 T applied along [001]. The inset 
shows the relative length change Al/l vs. H at T = 1.5 K 
and T = 0.2 K. The displayed curves represent raw data with- 
out background subtraction. The arrows indicate anomalies 
related to H C 2- 

FFLO state can exist in a type-II superconductor. They 
also calculated the dependence of the reduced tempera- 
ture to of a tricritical point on the Maki parameter a. 
For a = oo (orbital pair breaking is ignored), to = 0.55. 
Therefore, both Maki and FFLO scenarios give the same 
instability point for the second order phase transition in 
this limit. 

We do not see evidence for a second phase transition, 
which should be first order between FFLO state and a 
vortex state in FFLO scenario. However, we will make 
an assumption that we can use the result of the GG cal- 
culation of the instability point of the second order pase 
transition for arbitrary a within FFLO picture, to repre- 
sent the instability of the second order phase transition 
within the Maki's scenario for arbitrary a, since both 
theories give the same result in the limit a = oo, where 
calculations for Maki's scenario exist. To estimate a for 
CeCoIns we use Ao = 2.14fcsT c for a d-wave supercon- 
ductor and obtain H p = 2.25T C T/K = 5. 21. The orbital 
critical field is H c20 = 0.7H' c2 T c = 13.2 tM We there- 
fore obtain a — 3.6, corresponding to to = 0.35 from 
Fig. 1, curve b of Ref.Q. Alternatively, it is possible to 
find Hp from curve a of Fig. 2 of Ref. [| which relates 
H c2 (T = 0)/H p to a = V2H c20 /H p . Experimentally 
determined H c2 — 4.95 T, and the only unknown param- 
eter is Hp. We find H p = 5.8 T, a = 3.2 and t = 0.33. 
Both values of t are very close to the value of t = 0.31 
observed experimentally. This agreement leads us to con- 
clude that the first order phase transition is indeed due to 
the Pauli paramagnetism, confirming Maki's prediction 
of almost 40 years agctm. 

The FFLO state has also attracted great attention, 
but its unambiguous observation has not been made. In 
the last decade the FFLO state-was suggested to epst 
in heavy fcrmion UPd 2 Al 3 (Ref la) and CeRu 2 (RefM), 
based on thermal expansion and magnetization data, re- 
spectively. Subsequent research identified the magneti- 



zation feature in CeRu2 as due to the flux motionc3, and 
the region of the suggested FFLO state in UPd&Al 3 was 
shown to be inconsistent with theoretical modelEil. Most 
notably, no indication of the FFLO state was ever ob- 
served via specific heat measurement, which is the pri- 
mary tool for identification of the thermodynamic details 
of the phase transition. CeCoIns, according to theory, 
has all the prerequisite properties for the observation of 
the FFLO state, including that superconductor must be 
in clean limit. CeCoIns satisfies this requirement as well, 
with quasiparticle mean free path l tr ~ 14£o at T c t3. 
Within superconducting state thermal conductivity di- 
vided by temperature grows by an order of magnitude 
as temperature is lowered to T = 0.2T C , and even at 
30 mK, (1% of T c ) CeCoIns is outside of the impurity 
bandta. The agreement between theoretical prediction 
and experimental observation of the critical point at to 
shows that the pertinent physics of electronic spin Zee- 
man energy and Pauli effect drives the behavior of the 
system. However, within FFLO picture two transitions 
are expected: from normal state into FFLO state, likely 
to be second order, and from FFLO state into a usual 
Abrikosov vortex state via first order phase transition. 
We only observe one phase transition in CeCoIns with 
field H || [001]. Therefore, FFLO state is not present 
in CeCoIns, suggesting that FFLO state, if it exists, is 
much less stable than previously thought. 

In summary, the superconducting transition in 
CeCoIns with field H || [001] becomes first order below 
To = 0.31T C . This is consistent with long standing the- 
oretical predictions which take into account both orbital 
and spin interactions of superconducting electrons with 
magnetic field. We do not observe inhomogeneous super- 
conducting -FFLO state proposed theoretically almost 40 
years agooEl. Instead, first order superconducting tran- 
sition in CeCoIns takes place from normal metal into 
mixedryortex state, in accord with scenario suggested by 
MakBtJ. 
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